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Abstract. We introduce weighted succession rules and parametric production matrices — 
simple extensions of the standard ECO method succession rules and production matrices. 
The purpose is to enumerate combinatorial objects with respect to several variables. We 
consider one main example, from the theory of Dyck paths. The path statistics primarily 
considered are peak height, rise height, and certain subwalks. Many classical sequences 
(such as the Catalan, Motzkin, Narayana and Schroeder numbers) are incorporated in this 
example. 



1. Introduction 

Recently Deutsch, Ferrari and Rinaldi, [2], introduced production matrices as a device for 
ECO method calculations. Building on the concept of generating trees, the Enumeration of 
Combinatorial Objects method was developed and systematised by Barcucci and collaborators, 
[1]. It is a general method for enumerating combinatorial objects by studying how the objects 
grow according to a parameter. At the centre of method are the succession rules that describe 
how the objects grow — production matrices translate these to matrix notation. 

The production matrices considered by Deutsch et al. are integer valued. The step to in- 
clude matrices with parameters is natural and straightforward, and corresponds to weighted 
succession rules. In this paper we study a rich example of a parametric production matrix. 
Although only one matrix, it embodies many well known combinatorial number sequences, 
such as the Catalan, Motzkin, and Narayana numbers. The underlying objects for our ma- 
trix are Dyck paths, and the main statistics considered are heights of peaks and rises, and 
occurrences of certain subwalks. We recover known results about the height distribution of 
peaks, and use the matrix to show that it is equivalent to a rise height distribution. 

In Section 2 we give a brief description of the ECO method, define succession rules and 
production matrices, and introduce weighted successions. We recall results from [2] for op- 
erations on product matrices (which continue to hold in the parametric case). In Section 3 
we define Dyck paths and some of their statistics. We then define a parametric production 
matrix and show how it is related to those statistics. 

2. The ECO method 

The core of the ECO method is a recursive description of a class of combinatorial objects. 
This should be done in such a way that, if On denotes the set of objects of size n, each object 
O' G On+i is achieved from one and only one object O G On- We say that O' is a successor 
of O. 
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We assign a label (k), k G N"^ (the positive integers), to each object. An object's label 
give the number of successors of that object. The succession rule dictates the labels of these 
successors. The rule also includes an axiom (a), a G N"*", which gives the label of the smallest 
object. In the basic case a succession rule O is written as 



where ei{k) : N"*" N"*" gives the labels of the k successors of an object with label k. 

It is natural to think of a succession rule as a generating tree: the root of the tree is the 
axiom, (a), and if a node has label (k), it has k children labelled ei{k), 62 (fc), . . . , ek{k). 

A succession rule defines a sequence of positive integers an = the number of nodes at 
level n, n > 0. The generating function for a„ is defined by fn{z) = dnz""- 

The production matrix P{^) = P = {Pi,j)i,j>i of a succession rule is defined as follows. 
List all labels of a rule. Let this list be /i, /2, • • •■ Then equals the number of successors 
with label Ij produced by an object with label U. 

Example 1. The succession rule 



since label (1) has one successor with label (2) and (1) has two successors with labels (1) and 
(2). The sequence defined by Q. is the Fibonacci sequence (AOOOO45I), an = 1,1,2,3,5,8,.... 

2.1. Parametric production matrices and weighted succession rules. Naturally, one 
would like to use the ECO method to obtain not only the size distribution of the objects, but 
also multi-dimensional distributions. 

One strategy is to choose the labels wisely, since the distribution of labels on objects of 
a given size might be acquired. Indeed, using production matrices this is essentially as easy 
(or hard) as getting the size distribution. Another strategy is to use the succession rule and 
other data to find functional relations for the multivariate generating function. 

One of the aims of this paper is to demonstrate an alternative way of using production 
matrices to find multivariate distributions, by introducing weighted succession rules and para- 
metric production matrices. The idea is simple: 

Assume we are interested in the distribution of property, jointly with size. If a succession 
{k) (l) adds one property, we give that succession weight x. The sequence generated by 
the succession rule now is the weighted sum over objects of given size, and the generating 
function for the sequence is the bivariate generating function for property and size. 

Let W be a set of weights, and consider a matrix P = {pij)ij>i where pij S W. This para- 
metric production matrix corresponds to a weighted succession rule. A weighted succession 




= < (1) 



(2) 

(1)(2) 



has production matrix 




Six-digit numbers prefixed by 'A' indicate the corresponding entry in The On-Line Encyclopedia of Integer 
Sequences [4]. 
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rule is written as 

{a;w{a)) 



(k) {ei{k)-w (ei(fc))) (e2(fc); w (ei(fc))) • • • {ek{k)-w (ei(fc))) 



where again ei{k) : N"*" gives the labels of the k successors of an object with label 

/c, and w {ei{k)) G W gives the weight of that particular succession. In terms of generating 
trees, each edge is given a weight from W. We also have the possibility to give the axiom a 
a weight w{a). Of course, if W = {1} we recover the unweighted succession rules. 

From a weighted succession rule il, a parametric production matrix -P(O) = P = {pi,j)i,j>i 
is defined as follows. List all labels of a rule. Let this list be li,l2, ■ ■ ■■ Then pij equals the 
sum of weights of successors with label Ij produced by an object with label /j. The sequence 
{(J'n}n>o produced by 0, is now a sequence of polynomials in W, and we write 

oo 

/o(>V,z) = J^a„z". 

n=0 

Example 2. The succession rule 

0=<^(1)--(1;1) 

^(3)-(l;x)(2;l) 

has production matrix 

The sequence defined by $7 gives the Fibonacci polynomials (A011973) a„ = 1, 1, 1 + x, 1 + 
2x, 1 + 3x + x^, 1 + 4x + 3x^, . . .. 

Let us fix some notation. Let vJ = {1 . . .}, e = {1 1 1 . . .}"'' and / define a row vector, 
column vector, and identity matrix, respectively, of appropriate size. If a succession rule 
defines a sequence a„ with generating function /n(W',2;), we say that P = P{0,) defines a„ 
and write P ^ an and P — > /r2(W, z). 

We now cite some useful results for production matrices. These are from [2], where they 
were stated in terms of integer valued production matrices. To verify that these still hold in 
the parametric case is trivial. 

Fact 1 (Fact (iii), [2]). If P ^ an then an = P'^e. 

Fact 2 (Fact (v), [2]). If P ^ /(W,z) then f{W,z) = {I - zPy^e. 

Theorem 3 (Theorem 3.2, [2]). Ifh, c and r are nonnegative integers or weights and P — > 
f{W,z),then 

^_ , b ru^ \ ^ l + rzf{W,z) 



ce P J I - bz - rcz'^f{}V, z)' 
Corollary 4 (Corollary 3.1, [2]). Ifb, c and r are nonnegative integers or weights and 

then fiyV^z) satisfies 

rcz'^f^W, z)-{l-{b + r)z)f{W, z) + l = 0. 
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3. The Main Example 



3.1. Dyck paths. A Dyck path is a lattice path in the first quadrant from (0,0) to (2n,0) 
that uses steps (1, 1), we cah these up, or u, steps, and (1, —1) (down, or d, steps). We will 
call n the semilength of the path. Let D„ denote the set of Dyck paths of semilength n. The 
number of Dyck paths of semilength n are given by the nth Catalan number Cn (A000108), 



A peak in a Dyck path p is an occurrence of an up step immediately followed by a down 
step. The height of the peak is the height of the start point of the u step. A high peak is a 
peak at height 1 or higher. A rise is a maximal sequence of consecutive u steps. The height 
of the rise is the number of u steps minus 1. A high rise is a rise with positive height. A 
contact is any place the path touches the x-axis. A segment is any subpath consisting of 
consecutive steps. The segment is a level segment if it starts and ends at the same height. An 
excursion is the segment of the walk between two contacts with no other contact in between. 
Note that a Dyck path p with k + 1 contacts can be decomposed uniquely into its k excursions 

P = 6162 • • • efc. 




Figure 1. A Dyck path of semilength 9. The height of the rises are, in order, 
1,2,0,1,0. The path has 4 high peaks, of heights 1,3,3,3. Using succession rule 
of Theorem 5 the path will have weight a^o^^ix^yiylyi. Using succession rule 
of Theorem 8 the path will have weight yly^- 

3.2. The matrix. We are interested in the following matrix. 




The generating function for the Catalan numbers satisfies 



C{z) = l + zC\z)=Y,CnZ'' = - 



2z 



n=0 




Xo 
Xiy2 
X2y2 
X3y2 

xm 

xm 










Xo 

xiVi 

X2yA 

xzVi 







Xq 
X2y5 








Xo 

xiye 



\ 



p = 



X2yi 
x-m 
xm 

xm 



Xo 

xiVs 

X2y3 

xsVs 
xm 



V 



It is, as is shown below, the production matrix for Dyck paths, where the yk marks the 
height of rise number k and xj. marks level segments of the form u{uqkd)^d, where Qk are 



PARAMETRIC PRODUCTION MATRICES AND WEIGHTED SUCCESSION RULES 



5 



Dyck paths. We note that many classical number sequences arise from specialisations of the 
sequence associated with P. See Section 3.3. 

Theorem 5. For a Dyck path p, let hk{p) denote the height of rise number k in p and let 
Sk{p) denote the number of level segments of type u{uqkd)^d, where qk are Dyck paths. Then 
the weighted sum a„ over Dyck paths of semilength n is given by 

a„ := x^,°^\l'^''^x^,'^''\'^'^^ ■ ■ ■ = u'^P^e, 
and the generating function is given by 




(2) (2)x,y, (3) x„ 



Figure 2. The generating tree associated with the succession rule of Theorem 5. 

Proof. We will describe a weighted succession rule for Dyck paths whose parametric pro- 
duction matrix is P. The labels are the number of contacts (= the number of excursions + 
1). The axiom is the path ud with label (2;xo). The generating tree for this rule is sketched 
in Figure 2. 

Given a Dyck path p of length n with k excursions, p = 6162 • • • we now explain how the 
successors are derived. There are /c + 1 successors to p, with 1, 2, • • • , A; + 1 excursions. They 
are, respectively, 

1: ueie2 ■ ■ ■ Ckd 
2: eiue2 ■ ■ ■ Ckd 

k: 6162 • • • Ck-iuekd 
k+l: 6162 • • • Ckud. 
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Now, let a succession from a path with k excursions to one with I, 1 < I < k + 1, have weight 
yiXk-i+i- The succession increases the height of rise number / by one and adds one excursion, 
by default a level segment, of type u{uqdY^^^^d to the path as desired. Furthermore, each 
path has exactly one parent: the path obtained by by lowering the last excursion uqd to a 
Dyck path. □ 

In the following, we will use the notation x = (xq, xi, X2, . . .) and y = {yi,y2,y3, ■ ■ ■)■ If 
Xfc = 1 for k > n we will write F[xi, . . . , Xn', y, z ) instead of F{x]y; z), and similarly for y. 

Corollary 6. 



Fo .- F(xo, XI, , z) 

Proof. It follows from Corollary 4 that 

XQZ^Fi - (1 - (xo + xi)z) Fo + 1 = 0. 



Applying Theorem 3 we get 
Corollary 7. 

1 + xqzFq{xo,xi] ; z) 



□ 



Fi := F(xo,xi;yi;z) 



1 - xiyiz - XQyiz'^Fo{xo, xi; ; z) ' 
and in general 

1 + XQzFk_i{xQ,xi] ?/2, . . . , yt; z) 



Fk := F{xo,xi;yi,y2,...,yk;z) 



1 - xiyiz - xoyiz'^Fk-iixo,xi;y2, ...,yk;z)' 
Note the use of Fk-i{xo, Xi;y2, ■ ■ ■ ,yk] z) and not of Fi{xo,xi;yi, . . . ,yk-i;z). 

Below we find bivariate generating functions F{1, . . . , l,Xk', l',z) for length and height of 
rise number A: in a Dyck path. First we show that the heights of high rises have the same 
distribution as the number of high peaks. The reason for not using the perhaps more familiar 
high peak statistics in the definition of the main production matrix is that the joint distribu- 
tion of number of high peaks and level segments is not the same as the joint distribution of 
heights of rises and level segments. 

Theorem 8. The number of Dyck paths with whose high rises have heights (in order) ki, . . . ,ki 
equals the number of Dyck paths with ki, . . . ,ki high peaks of heights 1, . . . ,1, respectively. 



Proof. The result can be proved by describing an appropriate succession rule fi'. This will be 
such that -P(r2') = P{il.)\x=o. The generating tree for this rule is sketched in Figure 3. 

Paths are labeled by the number of steps in the last descent minus 1, except if that number 
is - then the label is the number of successive uds at the end, which equals the parent's 
label plus 1. The axiom is the path ud, and is given weight 1. The successors of a path p with 
label (k) are paths pi, . . . ,Pk+i with labels (1) to (/c + 1). The path pk+i is pud. The path pi, 
I < k, is achieved from p as follows. It is a combination of the paths p and q = u'+^d'"'"^; the 
latter is placed such that it ends at 2n + 2 and in pi we follow p up to the first intersection 
of p and q, after which we follow q. See Figure 4 for illustration. 
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(2) (2) y, (3) 1 



Figure 3. The generating tree associated with the succession rule of Theorem 8. 

Now, since a path's label gives either the number of uds at the end of p, or 1 plus the 
number of d steps in the last descent, the high peaks of p are preserved in pi, and if, I < k, pi 
has an extra peak at height /. 

It is thus clear that if we give the updates of (k) weights yi, . . . ,yk,l the succession rule 
corresponds to heights of high peaks as desired. It is also easy to see that each path has a 
unique parent. □ 




Figure 4. Example of the use of the concatenation rule used in the proof of 
Theorem 8. 

The next result is essentially equivalent to Theorem l.l(ii), [3]. Mansour gives generating 
functions for k peaks of height h in terms of Chebyshev polynomials. These can be summed 
to give an equivalent function for peaks of height h. As a function of t, the function G„ below 
is of the form and so can easily be expanded in t to give one-dimensional generating 

functions if desired. 

Corollary 9. Let Gn{t, z) be the generating function for high peaks of height n > 1, 

oo oo 

Gn{t,z) = 

i=0 j=0 
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where gij is the number of Dyck paths of semilength i + 1 with j high peaks at height n. 
Then' 

l-tz-tz^C{zr 

and for n >2 

z)zU{n - 4, - U{n - 2, 



Gn{t,z) 



Pn 



Gi{t, z)z^U{n - 2, - zU{n, ^) 

Proof. The first statement follows from Corollaries 6 and 7. 

Let Pn be the matrix P with all parameters except Hn set to 1. Now 

1 

e Pn-l 

By Theorem 3 the generating function for Pn^n > 2 is 

^ 1 + zGn-i 

1- Z - Z^Gn-1 

GizU{n-A,^)-U{n-2,-^) 

Now, straightforward but cumbersome calculations show that — :; — ^777 — , — — sat- 

Gi22c/(n-2,2-^)-C/{n,2^) 

isfies the recursion. □ 

Remark 1. The stated formula for Gn in fact holds for n = and n = 1 as well. 

3.3. Examples. As mentioned earlier, many well-known combinatorial sequences are con- 
tained in the sequence produced by the main matrix P. Here we list a selection of the most 
interesting of these. 

• F{] ; z) — > 1, 2, 5, 14, . . . = the Catalan numbers = the number of Dyck paths. (A000108). 

• -F(l, 0; ; z) — > 1, 1, 2, 4, 9, . . . = the Motzkin numbers = the number of Motzkin paths. 
(A001006). 

• F{2; ; z) — > 1, 3, 11, 45, . . . = the little Schroeder numbers = the number of Schroeder 
paths with no peaks at height 1 (A001003). 

• ; z) 1, 1 + t, 1 + 3i + t^, 1 + 6t + + t^, . . . = the Narayana numbers = the 
number of Dyck paths with k peaks. (A001263). 

• F(2s; ; z) 1, 1 + 2i, 1 + 6t + At^ , 1 + 12t + 24*^ + 8*^, . . . = the number of double 
rise bi-colored Dyck paths with k double rises. (Al 14687). 

• F(l + s; ; z) ^ 1, 2 + i, 5 + 5i + 14 + 2lt + + ^^ . . . = the number of Schroeder 
paths with k peaks, but none at height 1 = the number of double rise bi-colored Dyck 
paths with k double rises of a given color. (A126216). 

• F(l,2;;z) — > 1,3,10,36,137,... = the number of Schroeder paths with peaks at 
odd level. (A002212). 

• F(l, s; ; z) 1,1 + t,2 + 2t + t^ ,A + Qt + ?,t^ +t^ , . . . = the number of Dyck paths 
with k peaks at even height. (A091869). 

• F{t, 2; ; z) ^ 1, 2 i, 4 5i -h ^^ 8 18t + + t^,... = the number of hex tree by 
edges and left edges (A126182). 

• F(l, l + t]]z) ^l,2 + t,f> + At + t^,l'\. + Ibt + Qt^ + t^,... = the number of hex tree 
by edges and median children (A126181). 

• F(; 2; z) 1,3, 10, 35, . . . = the number of times a fixed Dyck word of length 2k occurs 
in all Dyck words of length 2n + 2k= i^^^^l), (A001700). 
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• F{;1 + t; z) ^ 1,2 + t,5 + At + 14 + 14t + Gt^ + t^, . . . = the number of leaves at 
level k+1 in all ordered trees with n+1 edges (A039598). 

• F(0,l;l + t;z) l,l + t,l + 2t + t^,l + ?,t + 'it^ + 1,... = Pascal's triangle (A007318). 

• F(0,l + t;l;z) l,l + t,l + 2t + t^,l + 2,t + ^f + l,... = Pascal's triangle (A007318). 

• F(l, 0; t]z) 1, 1, 1 + t, 2 + 2i, 4 + 4t + . . . = the number of Motzkin paths by 
height of final descent (A098979). 

Note that in all these cases Corollary 4 can be use to find an quadratic equation for the 
generating function. Conversely, we have the following result. 

Corollary 10 (Corollary 3.2, [2]). If a generating function f satisfies an equation of the form 

xoViz^f"^ - (1 - {xiVi + xo)z) / + 1 = 0, 

then f = F{xo,xi;yi; z). 
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